Critical view of the claimed 0+ pentaquark. 
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Abstract. We use a theoretical model of the y d — > K^K^ n p reaction adapted to the experiment done at LEPS where a 
peak was observed and associated to the ©+(1540) pentaquark. The study shows that the method used in the experiment to 
associate momenta to the undetected proton and neutron, together with the chosen cuts, necessarily creates an artificial broad 
peak in the assumed A'+n invariant mass in the region of the claimed ©+(1540). It is shown that the LEPS fit to the data, 
used to make the claim of the ©+(1540), grossly distorts the background. An alternative fit, assuming a background plus 
a fluctuation, returns a background practically equal to the theoretical one and a fluctuation identical to the one seen in the 
experimental p spectrum of 2cT significance. 
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In Ref.[l] the 7 '^C — > K^K X reaction was studied 
and a peak was found in the K^n invariant mass spectrum 
around 1540 MeV, which was identified as a signal for a 
pentaquark of positive strangeness, the "0+". The unex- 
pected finding lead to a large number of poor statistics 
experiments where a positive signal was also found, but 
gradually an equally big number of large statistics exper- 
iments showed no evidence for such a peak. A compre- 
hensive review of these developments was done in [2], 
where one can see the relevant literature on the subject, 
as well as in the devoted section of the PDG [3]. 

More recently a new experiment was done at LEPS 
on a deuteron target, and with more statistics, and a 
clear peak was observed around 1526 MeV in the K^n 
invariant mass distribution [4]. Yet, the experiment of [5] 
dealing with the same reaction as in LEPS but with ten 
times more statistics and with complete kinematics (but 
excluding small angles), failed to see any peak around the 
"0+" region. The detail of complete kinematics should 
be stressed because in the LEPS experiment neither the 
proton nor the neutron were measured and an educated 
guess had to be made for their momenta. 

In order to understand what is behind the peak seen 
in [4], we have constructed a theoretical model which 
contains the basic ingredients seen in the LEPS exper- 
iment, (p production on the proton and the neutron and 
A(1520) production on the proton, together with rescat- 
tering of the kaons. The details of this model can be seen 
in [6] . We adapt the model to the set up of the LEPS ex- 
periment, generating twenty random energies between 2 
GeV to 2.4 GeV and implement the angular cuts and the 
mass cuts to eliminate the (j) peak. 

The LEPS detector is a forward magnetic spectrome- 
ter. Its geometry is implemented in our simulation by im- 



posing that the angle of the kaons in the final state with 
respect the incident photon is not bigger than 20 degrees. 

The nucleons are not detected at LEPS, therefore, 
some prescription is required in order to estimate the mo- 
mentum of the p and n in the reaction yd K^K^np 
and determine the invariant mass of p or K^n. This is 
done using the minimum momentum spectator approxi- 
mation (MMSA). For this purpose one defines the mag- 
nitude 

Ppn = Pmiss = Pr + Pd- Pk+ -Pk- ( 1 ) 

which corresponds to the four momentum of the outgo- 
ing pn pair From there one evaluates the nucleon mo- 
mentum in the frame of reference where the pn system 
is at rest, pcM- Boosting back this momentum to the lab- 
oratory frame, we will have a minimum modulus for the 
momentum of the spectator nucleon when the momen- 
tum pcM for this nucleon goes in the direction opposite 
to Pmiss- Thus, the minimum momentum, /?„„■„, is given 
by 

_ 1^ I Emiss „ IPm/iil 
Pmin — - \PCM I ■ — H -fccM • — \^) 



where E^m ~ y \pcm\'^ + Mj;^ is the energy of the nu- 
cleon in the CM frame. In this case, the momentum of 
the other nucleon will be in the direction of the missing 
momentum with a magnitude 

Pres \Pmiss\ Pmin (3) 

In [4] the Mjf+„ invariant mass for the reaction yd ^■ 
K^K^np is evaluated assuming the proton to have a mo- 
mentum p,nin (actually what one is evaluating is Mj^+f^, 
with the non spectator nucleon). Consequently, in this 



prescription, the momentum of the neutron in the final 
state will be 

Pmiss 



Pn 



(4) 



which is used to calculate the Mr-+„ invariant mass for 
the reaction yd K^K^np in [4]. A cut is imposed at 
LEPS demanding that \pmin \ < 100 MeV. This condition 
is also implemented in our simulation of the process. 

In order to remove the contribution from the (j) pro- 
duction at LEPS one considers events which satisfy that 
the invariant mass of the K^K pair is bigger than 1030 
MeV and bigger than the value obtained from the follow- 



mg expression 



1020 MeV 



^0.09 X (£y^'^(MeV)- 



2000 MeV) (5) 



where Et^^ is defined as the effective photon energy 



2M„ 



(6) 



with Si^+i^-^ the square of the total center of mass energy 
for the K^K^n system calculated using the MMSA ap- 
proximation to determine the momentum of the neutron 
assuming the proton as spectator In [4] only events for 
which 2000 MeV < E'/^ < 2500 MeV are considered, 
a condition which is also incorporated in our simulation. 
The Ey^^ of Eq. (6) with the MMSA prescription would 
correspond to the photon energy in the frame where the 
original non spectator (participant) nucleon is at rest. 

The association of the K^N spectrum to K^n at LEPS 
has a repercussion in the assumed experimental K^n 
distribution. The real distribution from the model has no 
peak of the distribution around the "0+" peak. Instead 
we show in Fig. 1 the distribution obtained using the 
LEPS cuts and the MMSA prescription, normalized to 
the experimental data, which are also shown in the figure. 

We use the full model of [6] including production on 
the proton and the neutron plus A(1520) production on 
the proton. The MMSA prescription is used to compare 
to the data as explained before. As we can see in Fig. 1 
(histogram), the combination of the LEPS cuts and the 
MMSA prescription, which also affects the cut, has pro- 
duced an artificial peak below the region of the "0+". 
Coming back to the comparison of the "exact" distribu- 
tion, generated from the theoretical model, with the data 
(taken from Fig. 12 a) from [4]), we see that the "0+" 
peak has three points on top of the "exact" distribution. 
The accumulated strength of these three bins over the 
"exact" curve is about 35 events. The question now re- 
mains: could this peak, measured from the "exact" curve, 
be a statistical fluctuation due to the limited number of 
events of the experiment (around 2000 events)?. A hint 
to answer this question is provided by the LEPS experi- 
ment in fig. 10 a) of [4], where a peak followed by a dip 
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FIGURE 1. M^+„ invariant mass distribution obtained with 
the MMSA prescription and same cuts than LEPS normalized 
to the data of [4]. The thick black sohd line represents the 
fit of LEPS. The grey sohd fine (blue online) is the fit with 
the fluctuation. The dotted fine represents the background of 
the LEPS fit, while the dashed line (following the histogram) 
corresponds to the background of the fit with the fluctuation. 



is seen on top of the assumed background in the p 
mass distribution. In order to compare with these data, 
we show in Fig. 2 the mass distribution for/T^/? obtained 
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FIGURE 2. 

the MMSA prescription and same cuts than LEPS normalized 
to the data of [4] (shown as dots). 

with our model and the MMSA prescription (actually the 
K^N distribution, as discussed above), together with the 
data (Fig. 10 a) of [4]). The agreement with experiment is 
good, like for the K^n distribution shown in Fig. 1, indi- 
cating that we have indeed a realistic model. The discrep- 



ancies in the large momentum tail can be cured adding a 
small contribution of the broad (r= 300 MeV) A(1800) 
production to account for the tail of the K~ p distribution, 
which also reduces a bit the peak of the A(1520) upon 
normalization to the data, and we found minor changes 
in the A'+ n distribution, less than 5 % in the region of the 
"0+" peak. Yet, the point we want to make is that there is 
a peak in the data around 1650 MeV with four points over 
the "exact" curve with an accumulated strength of about 
33 events, followed by a similar dip, a typical structure 
of a statistical fluctuation. Since for this distribution there 
are no discrepancies in the interpretation with respect to 
[4] we can take the background from there (solid line 
of Fig. 10 a) of [4]) and we also find about 40 events. 
This peak was not associated to any resonance in [4]. In- 
stead it was dismissed as a statistical fluctuation. In fact, 
the peak disappeared when a slightly different cut was 
made (see Fig. 16 (Right) of [4]). The interpretation of 
this peak and dip as a fluctuation, as done in [4], is fair, as 
this peak has about 2a significance over the background 
((7 is the statistical error of the data for which \/Nevents 
is taken in [4]). But the peak of the "©+" over the "ex- 
act" distribution has also about 35 accumulated events 
and consequently could also be considered as a statisti- 
cal fluctuation. 

Here we would Uke to make a more quantitative study 
of the statistical significance of the "0+" peak, comple- 
menting the discussion done in [7]. In [4] a best fit to 
the data was done assuming a background and a Gaus- 
sian peak in the "0+" region. The best fit with these 
assumptions provided a background of about 22 events 
per bin below the "0+" peak, as can be seen in Fig. 12 
a) of [4] and reproduced in Fig. 1. With respect to this 
background, the "0+" peak has a strength of about 5 a. 
According to this, the statistical significance of the peak 
would rule out the possibiUty of it being a statistical fluc- 
tuation. 

Conversely, after the theoretical evaluation of the 
background, our argumentation goes as follows: The ac- 
tual background below the "0+" peak is bigger than the 
one provided by the LEPS best fit, around 36 events per 
bin instead of the 22 assumed in [4]. This makes the 
strength of the peak with respect to the background much 
smaller than in the LEPS best fit. It also makes a larger 
and the statistical significance is now of about 2a, some- 
thing acceptable as a fluctuation, as in the case of the ex- 
perimental K^p spectrum. This can be made more quati- 
tative by making the fit shown in Fig. 1 by means of a 
background and a fluctuation, with a gaussian peak fol- 
lowed by a gaussian dip, as it corresponds to the fluctu- 
ation seen in the p spectrum in Fig. 2. This fit to the 
data has about the same than the LEPS fit. Remark- 
ably, the fit with this shape returns a background practi- 
cally identical to the one evaluated with the theoretical 
model. 



This argumentation about the significance of the peaks 
is corroborated by further calculations. Indeed, in order 
to know the actual errors of the limited LEPS statisti- 
cal we have made 10 runs with the Von Neumann re- 
jection method, producing about 2000 events each, like 
in [4]. This method proceeds hke the experiment, gen- 
erating events or not according to their probability to be 
produced, and we have checked that the statistical signif- 
icance of the runs is equivalent to that of the experiment. 
From these runs we have evaluated the statistical errors 
of each run for each bin (see section VI of [6]). The er- 
rors found are of the order of 20% in the region of the 
Tjjis means an error (a) of about 7 events per bin, 
such that the difference of the peak to the background is 
indeed of the order of 2a. 

The other point we want to make is that it is possible 
to produce peaks by changing the cuts. This is shown 
explicitly in the experiment of [4] since the mentioned 
peak at 1650 MeV in the p mass distribution in Fig. 
10 a) of that paper disappears when a slightly different 
cut is made in Fig. 16 (Right) of the same paper. 

In summary, our study has shown that the background 
in the yd K^K n p reaction is fairly larger than the 
one obtained in the best fit to the data of LEPS assuming 
a background and a Gaussian peak in the region of the 
"0+". We also mentioned fliat flie fit of LEPS is not 
unique and other fits to the data, assuming a background 
and a fluctuation, are possible, producing the same re- 
duced and returning a background nearly identical to 
the calculated one. Based on the calculated background 
and the errors obtained from different Monte Carlo Von 
Neumann runs, we evaluated the statistical significance 
of the "0+" peak and found it to be of about 2a with 
respect to the background, compatible with a fluctuation. 
The larger statistical significance claimed in [4] was tied 
to the assumption of a significantly smaller background, 
which we have found is not justified. 
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